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PERTURBATION OF A GLOBALLY STABLE EQUILIBRIUM: 
APPLICATION ON AN AGE-STRUCTURED MODEL * 

Yannick Tchaptchie Kouakep** 


Abstract 

In this manuscript we consider an age structured epidemic system 
modelling the dynamics of transmission of immunizing disease like Hep¬ 
atitis B virus. Our model takes into account age as well as two classes of 
infected individuals (chronic carriers and acute infected human). Based 
on the low infectivity of chronic carriers, we study the asymptotic be¬ 
haviour of the system and, under some suitable assumptions, we prove 
the global stability of the endemic equilibrium point using perturbation 
arguments of semiflow. The intuitive conservation of the global stabil¬ 
ity under perturbation, needs in fact very technical tools and complex 
results. 

AMS subject classification: 35Q92, 34K20, 92D30. 

Keywords: Age-structured model, perturbation arguments, asymp¬ 
totic behaviour. Hepatitis B 


1. Introduction 

We consider two age-structured systems modelling an immunizing disease 
like hepatitis B and link them to the main partially ” aggregated” model to be 
studied in this manuscript. The first one uses the variables: 

a- s{t, a) the density of susceptible(s) at time t with chronological age a; 

b- i{t,T) the density of infective(s) that will develop acute disease at time 
t contaminated since a time r ; 

c- j{t,T) the density of infective(s) that will not develop acute disease 
(asymptomatic carrier(s)) at time t contaminated since a time t (some¬ 
times referred to as “age of infection” which plays a significant role in 
infective compartments i and j) 
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and reads as 

{dt + da)s{t, a) = a) — Ao(t)s(t, a), t > 0, a > 0, 

s{t,0) = A, 

(1 1) i9t+dr)i{t,T) = -{n + -fi)i{t,T), t>0, T>0, 

{dt +dr)j{t,T) = - {^l+'ye)j{t,T), t > 0, T > 0, 

pOO pOO 

i(i,0) = Ao(t) / pi{a)s{t,a)da, j{t,0) = Xo{t) / pj{a)s{t,a)da. 

Jo Jo 

Here A > 0 is some constant entering influx, p > 0 is the natural death 
rate, 7 ^ and 7 e are the additional death rates due to the disease. We set also 
(for (II.11) - (II.SI) ! the positive constants: i/j = {p + Ji) and z/j = (/r + jj). In 
addition pi € L“( 0 , 00 ) is a given function such that 0 < pi{a) < 1 a.e. 
while pj{a) = l—pi{a). Function pj represents the age-specific probability to 
become a chronic carrier when becoming infected at age a. Function pi denotes 
the probability to develop an acute infection when getting the infection at age 
a f fT4l[T5] b We refer to Edmunds et al.[^ for more explanation on the age- 
dependence susceptibility to the infection. 

The force of infection Ao(t) (for (ll.ll) - (ll.3D ') is expressed as 

poo 

(1.2) Ao(t) = / {Pi{T)i{t,T) + Pj{T)j{t,T))dT. 

Jo 

Here in Piir) and /3 j(t) denote the contact transmission rates between 

acute infected or chronic carriers with age of infection r with susceptible indi¬ 
viduals respectively. 

Finally this first model is supplemented together with some initial data 
s( 0 ,.) = so(.) e L+( 0 ,cx)) 

(1-3) 2 

z( 0 ,.) = fo(.), j( 0 , ■) = Jo(-) with (fo,jo) e [A+(0,oo)] . 

The second model with chronological age is 

[dt + da + p] s{t, a) = -Ai(t, a)s{t, a), 

[dt -f 9a -I- {pi + p + LJi)] i{t,a) = Xi{t,a)pi{a)s{t,a), 

[dt -1- 9o -I- iyj]j{t,a) = Xi{t,a)pj{a)s{t,a), 

[dt + da + p] r{t, a) = pii{t, a), 

posed for t > 0 and a > 0 , wherein t denotes time while a denotes the chrono¬ 
logical age of individuals, that is the time since birth. Here s{t, a) denotes 
the age-specihc density of susceptible, j{t,a) and i{t,a) denotes respectively 
the the age-specific density of chronic carriers and acute infected individuals 
(that can be symptomatic or asymptomatic) for Hepatitis B virus (HBV) while 
r{t, a) denotes the density of recovered individuals from HBV acute infection. 

Parameter p > 0 denotes the natural death rate, vj := {pj + p + ujj) > 0 
and z/j denote (for (ll.4D - (ll.7D ') the exit rates associated to each infected class 
i and j. The term vi gathers immunisation rate pi > 0 due to acute infection, 
natural death rate p and a possible additional death uji due to the infection; 
while i 2 j corresponds to death of chronic carriers, that is natural death rate 
and an additional death rate due to chronic infection and its consequences. 
Here we leave out possible recovery from chronic disease. The term Xi{t,a) 
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corresponds to the age-specific force of infection and follows the usual law of 
mass-action, that reads as 

nOO 

(1.5) Ai(t,a)= / [[3i{a,a')i{t,a') + I3j{a,a')j{t,a)]da . 

Jo 

Here in (lESl), I3i{a,a') and /3j{a,a') denote the contact transmission rates be¬ 
tween acute infected or chronic carriers of chronological age a' with susceptible 
of chronological age a respectively. 

This problem 11.41 is supplemented together with the boundary conditions: 

{ s(t, 0) = A, (constant influx) 

*(^i 0) = j(^i 0) = 0, (no vertical transmission), 

7’(t, 0) = 0, (no immunity at birth), 

and initial data 

(1.7) s(0,a) = so(a), i{0,a) = io{a), j{0,a) = jo{a), r{0,a) = ro{a). 

Note that the r—components of the system decouples from the other and have 
therefore no impact upon the long time behaviour of the system. It will be 
omitted in the sequel. Here recall that the boundary conditions i{t,0) = 
j{t,0) = 0 correspond to omitting vertical transmission (see [Tl 1161 [37l 1141 

m)- 

In the absence of the disease, the evolution of the population density follows 
the following simple (chronological) age-structured equation: 

(dt + da + fj) s(t,a) = 0, s(t,0) = A, 

where ^ > 0 denotes the natural death rate while A > 0 corresponds to the 
birth rate, that is assumed to be a constant influx. 

With 350 million people chronically infected, the Hepatitis B virus (HBV) 
pandemic is a global public health problem and constitutes one of the main 
cause of cirrhosis and hepatocellular carcinoma with a high rate of morbidity 
and mortality especially in South-East Asia and Sub-Saharan Africa struggling 
with a prevalence of HbsAg (Hepatitis B surface antigens) chronic carriers 
greater than 8% |24111711121132) . The class of HBV chronic carriers has a high 
risk of developing liver disease such as the one mentioned above. 

In highly endemic areas the epidemiology of HBV exhibits two main fea¬ 
tures: a low average age at infection and a high prevalence of chronic carriers. 
According to numerous works (see for instance [H (S) 11611181125) and the ref¬ 
erences therein) there is a rapid decline in the probability of developing the 
chronic carrier stage of the disease with respect to the age at infection. Accord¬ 
ing Edmunds et al [^ , approximatively 90% of children of less than 5 years old 
will become carriers if infected while about 10% of adults will become carriers 
if infected. Moreover this age-specific probability appears to be remarkably 
stable for a wide range of areas [6]. Hence the mathematical modelling of 
HBV infection requires to take into account this strong age-specific differential 
susceptibility of humans. 

In such highly endemic areas, these are two main routes for the infection of 
children: a vertical transmission from an infected mother to her infant, and an 
horizontal transmission through close contacts. 
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Following WHO [3 7) . perinatal transmission is a much less important con¬ 
tributor to the carrier pool in Africa (we also refer to Assumption 4.1 in [38] 
for similar consideration for the Chinese case or fniTG) for the Pakistanese 
case). Moreover Anfumbom et al. state that “A low proportion of HBeAg 
among HBsAg-positive pregnant women with known HIV status could suggest 
low perinatal transmission of HBV for an endemic country like in Cameroon” 
(see also fT4llT5]l. Hence, as in the work of Edmunds et al. [610, we shall 
focus, in this work, on horizontal transmission taking into account age-specific 
susceptibility. 

This differential susceptibility is a particularly important point for HBV 
infection. Let us recall that according to CDC^ (see also |6] and the refer¬ 
ences therein) about 90% of children will remain chronically infected with HBV 
while 90% of adults will develop acute infection and will completely recover 
from HBV infection. 

Hence the main feature of the above described models consists in the age- 
specific susceptibility dependence through function pj (a) [mill] (or pj{a)). 
Edmunds et al. in [6| justihes the fact that “the expected probability” pj (of 
developing the chronic carrier state given the age a at infection in an expo¬ 
nential function of age a) is decreasing. More specifically according Edmunds 
et al. in [6] function q = pj takes the form pj{a) = pi{a) = for some 

suitable parameter set. In order to take into account this age-specific suscepti¬ 
bility dependence we will use in this work a simplest prototypical shape curve 
of the form 

(1.8) pj{a) = Ke“™, 

for some k S [0,1] and r > 0. Note that this specific form will only be used 
to check some assumptions in order to apply our main perturbation result. 
Based on the data of [6] we estimate the above parameters k and r using the 
least squares method to obtain the following function (graphically closed to 
Edmunds’slU Fig. 1, page 199] function pi for r = 0.645 and s = 0.455): 

(1.9) pj(a) = p 2 (a) = 0.643exp (—0.156a), a > 0. 

The above models (fTH-lfOl) and pm -dTTTl) have been suggested in the work 
of Bonzi et al using ordinary differential equations and discrete age structure 
for susceptible population. We refer to Zou et al [38] and the references therein 
for other classes of age structured system modelling HBV transmission. 

Standard methodologies apply to provide the existence and uniqueness of 
mild solution for Systems (II. 1111 - (II. 1311 and (ll.4ll - (ll.7ll (see for instance Theo¬ 
rem 12.31 in section [21 also [^[^[30] and the references therein). 

The aim of this manuscript is to investigate the asymptotic behaviour 
of an “aggregated” system concerning infective compartments (with I(t) := 
i{t,x)dx and J{t) := j{t,x)dx) derived from systems (ll.lll - (ll.3p and 
dTU-dni). To do so let us recall that according to WHO [36], Edmunds et al 
[9] , Bonzi et al [3] , Fall et al |llj and Wilson et al [341135] , chronic carriers 
(most of time asymptomatic) have a low infectious rate. As a consequence in 


^Centers for Disease Control and Prevention, USA: www.cdc.gov 
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most part of this work we assume that 

(1.10) 0 < p, « A. 

In the extreme case when A = 0 systems (El or El re-write as the 
so-called reduced system (omitting the recovered class). 

We shall assume that the contacts between individuals are homogeneous 
among the different cohorts, so that functions A and A are constant 

A = A > 0 and A = A > 0- 

We will then study the general semiflow generating (s, /, J) as a perturbation 
on Pj around 0, of the reduced systems. 

In addition for each function p € L°“(0, oo) we set p* G (l^^(0, oo)) the dual 
form associated to p and defined by 


p [p\ = 


p{a)(p{a)da, \/(p G L^{0,oo). 


Our technical main “aggregated” system (I1.11I) - (I1.13I) is obtained from 
both models im-lfOl) and (ira-lfTTl) by setting I{t) := i{t,x)dx and 
J{t) ■= j{t, x)dx, such that (s, I, J) satisfies the following closed system of 
equations: 

' (A + da+ p) s{t, a) = -X{t)s{t, a), t > 0, a> 0, 
s{t,0) = A, 

r{t) = \{t)p* [s(t,.)] - u//(t), 

J'{t) = A(t)p}[s(t,.)] - ujJ(t), t > 0, 
wherein we have set 


( 1 . 11 ) 


(1.12) X{t)=PiI{t) + PjJ{t). 

Next by setting /q := io(x)dx and Jq := jo{x)dx one supplements 
(11.111) together with the initial data: 

(1.13) s(0,.) = so(.) e ^^^(0,00), J(0) = Jo>0, J(0) = Jo>0. 

We investigate the long time behaviour of (II.Ill) - (II. 131) when Pj is small using 
some knowledge about the dynamics of the auxiliary reduced systems with 
Pj = 0. The dynamics for the chronic carriers is fully known from the dy¬ 
namics of the reduced system because the (s, J)—variables decouple from the 
J—variable. 

In order to understand the dynamics of the above problem (I1.11I) - (I1.13I) with 
0 < /3j < < 1, we will make use of the abstract results of Magal in uni to derive 
a result about global asymptotic stability of the endemic equilibrium for small 
values of Pj. These abstract results of Magal [20] developed in the persistence 
context, the global perturbation of stable equilibrium theory initiated by Smith 
and Waltman |28j . 

The work is organized as follows. In Section 2, we deal with the well- 
posedness of the system, derive preliminary results that will be useful to study 
the long term behaviour of the model and we study the global asymptotic 
stability of the disease free equilibrium in the case when the basic reproduc¬ 
tion number satisfies i?o < 1 (see (12.211) for the definition of this threshold 
parameter). In Section 3 we state and prove our main perturbation result. 
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It is concerned with the long time dynamics of (II.1111 - (II.131) with Pj « \ 
small enough. It is stated under rather general assumptions without using the 
specific form of function pj while the proof is based on the perturbation argu¬ 
ments derived by Magal in |20j . Finally Section 4 is devoted to an application 
of this perturbation result to our prototypical function (11.81) and concluding 
remarks end the manuscript. 


2. Preliminary results 

The aim of this section is to investigate basic properties of System (11.111) - 
(11.131) . This analysis will be related to the following assumptions: 

Assumption 2.1. We assume that A > 0, p > 0, /3i > 0, Pj > 0, i/j > 0 
and vj > 0 are given parameters. 

Our second assumption is related to function pj. 

Assumption 2.2. We assume that pj G TJf’(0,oo) is a given measurable 
and bounded function such that 

0 < pi{a) < 1, for a.e. a > 0. 

In order to deal with (I1.11I) - (I1.13I) . let us introduce the Banach spaces 
X = L^{0, oo) X M X R X R, Xq = T^(0, oo) x {0} x R x R 
endowed with the usual product norm, as well as its positive cone A_|_ defined 

by 

A+ = T+(0, oo) X [0, oo) X [0, oo) x [0, oo) and Xo+ = Aq D A+. 
Consider the linear operator A : D{A) C X ^ X defined by 


D{A) = oo) X {0} X R2 and A 


Consider also the non-linear map of the class C °°, F : Ao+ —>■ X defined by 

/ -{Piai + Pjaj)p \ 

A 

{Piai + Pjaj)p*j[ip] 

\{Piai + Pjaj)p*j[tp\) 


( 


/-ip' - pip\ 

0 


-p{0) 

ai 


-vjai 

\ajj 


\ -vjaj ) 


F 


linear 

map ( 

1 

/ 

0 


ai 


\aj) 

V 


Now identifying {s{t, .),I{t),J{t)) together with u{t) = {s{t, .),0,1{t), J{t)) , 
Systen 

(2.14) 


System (I1.11I) - (I1.13D re-writes as the following abstract Cauchy problem 
u'{t) = Au{t) + F {u{t )), t > 0, 
u(0) = (so(.), 0, Jo,-^o)^ € Ao+. 


Together with these notations, our first results are collected in the following 
theorem: 


Theorem 2.3. Let Assumptions \2.1\ and \2.‘A be satisfied. Then System 
(I1.11I) - (I1.13I) generates a strongly continuous positive semiflow {U{t)}t>o on 
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This means that for each x € the continuous map u : t ^ u{f) := 

U{t)x defined from [0,c») into Xo+ is a weak solution of (12.141) . that is 

f u{s)ds € D{A), Vt > 0, 

Jo 

u{t) = X + A / u{s)ds + / F{u{s))ds, Vi > 0. 

Jo Jo 

The semiflow {C/(i)}t>o satisfies the following properties 

(i) ifx = (so, 0, Iq, Jq)^ S Xq+j then denoting U{t)x = {s{t ,.), 0, /(i), J(i))^, 
it satisfies the following Volterra integral formulation: 

jsoia-t)e-J>+^^A]d^ifa>t>0 


while functions t —^ I[t) and t —> J(i) are of class on [0,cx)) and 
satisfy 


(l{0)=lo, J(0) = Jo, X{t)=(3iI{t) + f3jJ{t), Vi>0, 
(2.16) < r{t) = X{t)p*j[s{t ,.)] - 

[ J'(<) = A(i)p}[s(i,.)] - ujJ{t), t > 0. 


(ii) It satisfies the following bounded-dissipative estimates for each x £ Aro+ 
and each i > 0; 

\\U{t)x\\x, <^{\-e-'^) + \\x\\xe-'\ 

wherein we have set v = min {/i, u/, uj}. 

(iii) The semiflow {t/(i)}t>o is asymptotically smooth on Xq- 


Proof. The proof of (i) is rather standard. Indeed it is easy to check that 
operator A satisfies the Hille-Yosida property. Then standard methodologies 
apply to provide the existence and uniqueness of mild solution for System 
(I1.11I) - (I1.13|) . (see for instance [2n[22l[30] and the references therein). 

The proof of (ii) is immediate from the integration of the equations. 

We will now focus (iii). In order to prove this result we will make use of 
results derive by Sell and You in |26j . More precisely, we will show that for 
each bounded set B C Aro+, the semiflow {U(t)}t>o is asymptotically compact 
on B. Let B C Yo+ be a given bounded set. Let us consider a sequence 
(^si, 0, Iq , Jg^ ^ C B of initial data and let us denote by 

{s^t ,.), 0,P(t), P(t)f = U(t) (s^o,0, /j)', 

the solution semiflow. Let {tj}j>o be a sequence tending to oo. We aim 
to show that the sequence {s^(tj, .), 0, I^(tj), J^(tj))^ is relatively compact in 
Xq. To that aim, let us first notice that due to estimate (ii) in Theorem 
IQ the sequence of functions J^(tj + .) and I^(tj + .) are uniformly bounded 
as well as their time derivatives. Using the Arzela Ascoli theorem, possibly 
along a subsequence, one may assume that (J^(tj + .),P(tj + .)) converges 
locally and uniformly to some bounded, continuous and positive functions t £ 
K —)• (J(t),I(t)). Then setting for each j > 0, Xj(t) = l3iP(t) + PjJ^(t) and 
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A : R —[0, cxd) defined by X{t) = + /3jJ{t), one obtains using Lebesgue 

convergence theorem that as j —>■ oo 



a — tj)dae 


- lo^ [f<.+\j{s)]ds 


< 



sl{a)dae —>■ 0 . 


while as j —oo 


g-_ g- /“„[M+A3(t3+s)]<is g- /“„[/i+A(s)]ds 


for a.e. a > 0. 


As a consequence one obtains that 

and 

lim Sj{tj ,.) = in L^(0,oo). 

j->-oo 

The result follows. q.e.d. 


From Theorem (53] one deduces using the results of Hale |13) (see also Smith 
and Thieme m Magal and Zhao |23) and the references therein), the follow¬ 
ing results: 

Proposition 2.4. Let Assumvtion \2. be satisfied. The semiflow {U{t)}t>o 
provided by Theorem \2.S\ has a non-empty compact global attractor A C Alo+. 
This means that A is compact, invariant and attracts all bounded set B C Xq^, 
in the sense that for each B C Alo+ bounded subset, one has d(U{t)B,A) —>■ 0 
as t —>■ oo where d{B,A) denotes the semi distance from B to A defined by 

d{B, A) = sup inf \\y - x\\x. 

Due to the above result the description of the long time behaviour of (11.1111 - 
(11.131) relies on the suitable descriptions of the global attractor. Now we need 
to derive some basic properties of the entire solutions of (I1.11I) - (I1.13I) . Before 
doing so, let us recall that the entire solutions are smooth, in the sense that 
s G (^““(R X [0, oo)) while (/, J) G C°°(R)^. With this in mind, we can state 
our first lemma. 


Lemma 2.5. Assume that the assumptions of Proposition \2f4\ are satisfied. 
Let t G R !-->■ (s(t, .), /(t), J(t))^ be a given entire solution of (I1.11I) - (I1.13I) . 
then the following property holds true: 

(2.17) s{t,a) < Ae~^°' := spia), V< G R, V a > 0. 

One furthermore has: 

(2.18) (3to G R flilih) + fljJito) = 0) ^ (/(t) = J{t) = 0, Vt G R) 


and, 

(2.19) 

{3to G R, Oo > 0 s(to,ao) 


Ae-/^“o) => 


/(t) = J(t)=0,VtGR, 
s{t, .) = sf(.), t G R. 
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Proof. In order to prove (12.171) we make use of the integral Volterra formu¬ 
lation of the solution that reads as follow: for each t G K and a > 0 

(2.20) s(t, a) = 

Hence since A > 0, (12.171) follows. 

Let us now prove (12.181) . With this in mind, denote by 77 := max(L'/,uj) 
and V := min {vj^vj). Then due to (11.111) one has: 

K't) [git) -P]< A'(t) < X{t) [g{t) -e], Vt G R, 

where we have set g{t) := [/3/p/(a) -|- /3jpj(a)] s{t, a)da. Then if there exists 

to G R such that A (to) = 0 then due to the above differential inequality one 
obtains A(t) =0. As a consequence, since /?/ > 0, on the one hand if /3j > 0 
then the result follows. On the other hand, if /3j = 0 then one gets from the 
above that I{t) = 0 and J becomes a uniformly bounded entire solution of 
J'(t) = —vjJ{t) for all t G R. Hence J(t) =0. 

Finally to complete the proof of (12.191) . Note that if there exist to gM. and 
oq > 0 with s{to,ao) = spiao) then using (I2.20|) . there exists ti G (to ~ ooito) 
such that A(ti) = 0. Thus due to (12.181) . A(t) = 0 while (12.201) ensures s{t ,.) = 
s_f(.). The result follows. q.e.d. 


In order to give some information on the long term behaviour of system 
(I1.11I) - (I1.1.‘1() . we introduce the threshold parameter Rq defined by 


( 2 . 21 ) 


Rq A 


Pi [e 




Pj [e 


Using this definition, straightforward computations lead us to the following 
result: 


Lemma 2.6. Let Assurtmtion \2.R2.^ be satisfted. Then the following holds 
true: 

(i) If Rq < 1, then System (Il.lip - (I1.13I) has a unique stationary state 
xp = (sf, 0, 0,0)^ G Xp where spia) = Ae“^“. 


(ii) If Rq > 1, then system (Il.ll|) - (jl.l3|) has two stationary states xp G Xp 
and Xp = (sb(-)) 0: Ep)'^ with 

sp(a) = sp(a)e~^‘^°', Ip = —pjlsp], Jp = —p}lsp], 

Vi vj 

and where Xp > Q is the unique solution of the equation 


1 = A E 



^ — {p+Xe). 


Together with this threshold definition, we are now able to study the global 
dynamics of (II.Ill) - (II. 131) when i?o < 1. Our result is stated in the next 
proposition: 


Proposition 2.7 (The case Rq < 1: disease extinction). Let Assumption 
lif.lHil.ill be satisfied. Assume that Rq < 1, then the global attractor A defined 
in Proposition\2.4\ satisfies A = {tb}- 
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Proof. Let us consider for k S {I,J} the quantities 
Tk = —■ Using these notations, note that one can re-write 


> 0 defined by 


k£{I,J} 

Next consider the map L : R defined by 


L 




Til + TjJ. 


Let X & Ahe given and let us denote by |u(t) = (s(t,.), 0, /(t), J(t))^| C 

A an entire solution of (11.111) with V (0) = x. Then straightforward computa¬ 
tions yield the following: for each t S R and r > 0: 

( 2 . 22 ) 


L[V{tPT)\-L[V{t-T)\ 





Recalling that Rq <1 and using (12.171) . one obtains that the map t ^ L [U(<)] 
is bounded and decreasing on R. 

Let {tn}n>o be a decreasing sequence tending to —oo as n —>• oo. Consider the 
sequence of maps {Vn = U(. -|- tn)}n>o- Up to a subsequence one may assume 

that Vn converges towards some function V = 0, 1, for the topology 

of Cioc (R, ^). Note that |u(t) | is also an entire solution of (I1.11I) - (I1.13I) . 

L J tgB 

Using (12.221) . one obtains that V satisfies: 


(2.23) 


^ r,p:.[s(t,.)]-i 

fee{7,J} 


+ pjJit) 


= 0 . 


Furthermore since {t„} is decreasing and the map t ^ L [U(t)] is decreasing, 
this leads us to 


(2.24) 


0 < L [V{t)] < L U(0) 


Vt e R. 


Now if i?o < 1, we infer from (12.231) that I{t) =0 and J{t) = 0. From (I2.24L 
one gets that I{t) = J{t) = 0 and since V is an entire solution of (I1.11I) - (I1.13L 
one obtains that V(t) = xp and the result follows. 

If we assume that i?o = 1, then, by using Lemma 1^31 (12.181) and (I2.19p . (12.231) 
implies that either s(t,.) = sp or I{t) = J{t) = 0. The second condition can 
be handled similarly to the argument for i?o < 1. In the first case, namely if 
s'(t,.) = sp, one gets from the s—equation that 

(/3,/(t)+/3jJ(t))sF(a)=0, 
so that + PjJ{t) = 0 and the result follows. 


q.e.d. 
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3. The case i?o > 1: a perturbation result 


In this section we will study global dynamical properties of (jl.lll) - (ll.l3l) in 
the framework of the biological assumption << 1 being a small parameter. 
For notational simplicity we replace /3j by e. Moreover we will explicitly write 
down the dependence of several quantities with respect to e > 0. For instance, 
recalling (j2.21|l . we write Rq = i?o[e]^ U{t) = Ue{t) the semiflow provided by 
Theorem 12.31 xe = x% the endemic equilibrium defined in Lemma 12.61 and so 
on. 

Before stating our main result, let us state our assumptions that will be 
discussed in the next section using the specific form of function pj in (11.81) : 


Assumption 3.1. We assume that the following properties hold true: 

(i) i?o[0] > 1 and, 

(ii) The semiflow {t/o(0}t>o satisfies: 

\^Uo{t)x = XE, Vx e M* := |(so,0,/o, Jo)^ e A+: /q > o| . 

(iii) The map A : > C defined by 

= {z G C : Re ( 2 ) > —p} 


A{X) = X + lEl3h*i 


se{-) 



does not have any root with Re(A) > 0. 


Then the following result holds true: 

Theorem 3.2. Let Assumvtion \2.ll \2.S\ and \3.1\ be satisfied. Then there 
exists (5 > 0 such that for each e G (0, i5), the semiflow {f7e(t)}t>o satisfies 

lim Ue{t)x = x%, 

t —¥00 

for each x = (so, 0, Iq, Jq)^ G Ao+ such that Iq + Jo > 0. 

The proof of this result relies on the application of Theorem 1.2 derived by 
Magal in |20j (see also Smith and Waltman |28j for an earlier result). In order 
to use the above mentioned result, let us check the required set of assumptions. 
They are collected in the following lemma: 

Lemma 3.3. Let Assumvtion [KT\ \2.2\ and \3.1\ (i) be satisfied. Consider the 
map p : [0,00) x ^0+ ^ [0, 00) defined by 

p(e, x) = pe{x) = fill + eJ, Vx = (s, 0, 1, J)^ G Ao+, 

and let us set for each e > 0; 

(3.25) M|f = {x G Xo+ : Pe{x) > 0} , dM^ = {x G Xo+ ■ Pe{x) = 0} 

Then the following holds true: 

(i) For each e > 0, Mq and DMq are both positively invariant under the 
semiflow {f7e(t)}t>o, 

(ii) For each e > 0, Lf^ has a global attractor Ae in A-|_ and the family 
{■^e}E>0 is upper semi-continuous at e = 0. 
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(iii) For each 6o > 0, there exists rj > 0 such that for each e € [0, (5o] and each 
X G Mq, one has 

liminf (JJe{t)x) > rj. 

t—¥00 

Before proving this result let us first complete the proof of Theorem 13.21 As 
already mentioned the proof of this result is a direct application of Theorem 
1.2 in Magal |20j . To that aim, let us introduce for each e > 0 the map 
Fg : Xo+ —?► X defined by = F^+eG where Fq : Xo+ —> X and G : Xo+ —?• X 
are defined by 


( 


/ -(iiaup \ 


( 


/ -ajp \ 

0 


A 


0 


0 

cti 


Piaip} [p] 

, G 

ai 


ajP*i W\ 

\ajj 


\Piaip*j [p]) 


\aj) 


\ajPj M/ 


Let us introduce the linear semiflow {L{t)}t>o C C{Xo) defined for each 
x G Xq by L(t)x = u(t) where u{t) is the mild solution of the following linear 
Cauchy problem: 

(3.26) ~ Au{t) + DFq (^x%) u{t), t > 0 and u(0) = x. 

We first claim that: 

Lemma 3.4. For each t>0 the spectral radius of L{t) denoted by r{L{t)) 
satisfies: 

r{L{t)) < 1. 

The proof of this lemma is postponed. 


Let us denote for any subset A C Xqj^ and map f : A ^ X: 

\\f{x) - f{y)\\ 


(3.27) 

Using this notation we claim that: 


||/||Lip. a:= sup 

x,yGA, x^y ||^ VW 


Lemma 3.5, For any given and fixed io > 0 one has: 

(3.28) ^Ihn \\Ue{to) - L{to)\\up,Bx„{.%,s)nXo+ = «■ 

Before proving this lemma, let us notice that Theorem 13.21 directly follows 
from the properties stated in Lemma 13.31 Lemma 13.41 and Lemma 13.51 

The rest of this section is devoted to the proof of Lemma 13.31 Lemma 13.41 
and Lemma 13.51 


3.1. Proof of Lemma 13.31 This section is devoted to the proof of Lemma 

EH 

Let us first notice that (i) is straightforward. 

In order to prove {ii) we will make use of Proposition 2.9 in m Recall that 
due to Proposition 12.41 for each e > 0 the semiflow {Ue{t)}t>o is bounded, 
dissipative, asymptotically smooth and has a global attractor denoted by Af, C 
Aro+ that attracts all bounded subsets. According to Proposition 2.9 in |20j . 












PERTURBATION OF A GLOBALLY STABLE EQUILIBRIUM 


13 


in order to prove that the family {^e} is upper semi-continuous at £ = 0 it is 
sufficient to show that 


(3.29) Ue{t)x —>• Uo{t)x, as £ —>■ 0'*', 

uniformly with respect to (t, x) on bounded sets. This is an application of 
Gronwall’s Lemma. Let r > 0 and k > 0 be given. Set B = Bxo (0: l~l Xo+ 
and note that due to Theorem 13.21 (ii) one has for all £ > 0 and t > 0: 

Ue{t)B C B with B = Bxo ^0, — -b fl Xo+. 

Denote by Aq : D{Ao) C D{A) —>■ D{A) the part of A in Xq = D{A), which is 
defined by 

(3.30) Aqx = Ax, \/x € D{Aq) = |a; S D{A) : Ax € D(A)| . 

Recall that Aq is densely defined and generates a C^-semigroup denoted by 
{TAo(t)}t>Q in Xq. Here and in the sequel, we will denote by i?(A;H) = 
{XI — A)~^ the resolvent of A for each value A £ p{A), the resolvent set of A. 
Let X £ B he given. Then for each £ > 0 and f > 0 one has: 

Ue:(t)x = TAo(t)x + lim [ TAo{t- s)XR{X;A)Fs{Ue{s))ds. 

^-^00 7 g 

We refer to Magal and Ruan [22) for details on the above constant variation 
formula. Hence, since A satisfies the Hille-Yosida property, there exists some 
constant M > 0 such that for all t £ [0, r]: 

||17e(t)x - Uo{t)x\\ < eM sup \\G{y)\\ + MIIFoIIlip. b [ “ t^o(s)a:||ds. 

yGB *^0 

Then by using Gronwall’s inequality, one completes the proof of (I3.29|) and {ii) 
follows. 

Next it remains to prove {in). To that aim we will prove uniform weak 
persistence and will convert later into uniform strong Dersistence |31) . We 
include parameter e into the state space. We consider the Banach space Y = 
X X R and define: 


Yo = XoX R, Y+=X+x R+ and ^0+ = Yo G Y+. 
Then we consider the following problem: 

f u'{t) = Au{t) + FQ{u{t)) -b eG {u{t)) , t > 0, 


(3.31) 


e'{t) = 0, 


supplemented together with initial data (m(0),£(0)) G Yq^. Let us consider 
the semiflow {H(t)}t>o on Lo+ generated by (13.311) . It is defined by 


V{t) 


Ue{t)x 


V 


£Yc 


o+- 


Let (5o > 0 be given. Let us also consider the complete metric space M := 
Xo+ X [0, Jq], the continuous function p : M —>• [0, oo) defined by 


(3.32) 


p{y) = I + eJ, y = (s, 0, 1, J, e) £ M, 










14 


Yannick TCHAPTCHIE KOUAKEP 


as well as the sets 

(3.33) Mo = {y e M : p{y) > 0} and dM^ = M \ Mq. 

Then the following result holds true: 


Lemma 3.6. Under the above assumptions and using the above notations 
the semiflow {V(t)}t>o is p weakly uniformly persistent, namely there exists 
77 > 0 such that 

limsup p{V{t)y) > 77 , Vy S Mq. 

t—¥OC) 


Proof. In order to prove this result, let us first notice that due to Theorem 
12.31 for each y £ Mq there exists a time ty > 0 such that 


V{t)y £ Mo n 




X [0,(5o] 


Vt > ty. 


As a consequence to prove Lemma 13.61 it is sufficient to show that for each 
r > 0 there exists 77 ^ > 0 such that 


limsup p{V{t)y) > 77 ^, V 77 £ Mo(r-), 

t—^OC 

wherein we have set Mo(r) = Mq fl [Bxq (0, r) x [0, Jo]] ■ 

Let r > 0 be given. Let us assume by contradiction that there exists a sequence 
{Vn = {xn,£n)}n>o C Mo(r) and a sequence {tn}n>o tending to 00 as n, —>■ 00 
such that: 

p{V{tn + t)yn )<——, Vtt > 0, t > 0. 

77+1 

Up to a subsequence, one may assume that eo G [0, Jo] and possibly up 
to extraction of an other subsequence one easily checked for the system that: 

(3.34) lim V (t„ + t)yn = 7/00 

n^oo 

uniformly with respect to t > 0 . 

Next, as a consequence of (13.3411 . one obtains that for each «: > 0 small enough 
there exists n = n^n) > 0 large enough such that for any k £ {/, J}: 

(3.35) pI[sf] - K ^plls'^it,.)] <pI[sf] +K, Vt>0. 



Here we have set for each t > 0 and n > 0. 

(s"(t,.), 0, /'‘(t), J^{t),enf = U(t„ + t)77„ 


Note now that for each n > 0 and each t > 0 one has: 


(3.36) (/3iM(t) + SnJ^it)) > ^iM(t). 

Let K > 0 be given such that p} [sf] — k > 0. Let n = n^n) be defined by 
(j3.351) . Then we infer from ()3.35|1 . ()3.36(1 and the /—equation that for the 
above defined value of n, namely n = n^n), one has for each t > 0 

(/”)'(t) > /3//”(t) [p*i [sf] - k] - ^'/-^"(f), 

This re-writes as 

/"(t) > /”( 0 )e^‘, Mt > 0 , 
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as soon as k > 0 is small enough, namely ^ (i?o[0] — 1) > R- Here recall that 
one has assumed i?o[0] > 1 in Assumption 13 . iH il . Since /”(0) > 0 this last 
sub-estimate contradicts the boundedness of I. This completes the proof of 
Lemma 13.61 q.e.d. 

Let us recall that the semiflow {H(t)}t>o is bounded dissipative on M and 
it is asymptotically smooth so that the semiflow has a compact attractor that 
attracts all bounded subset of M. (see Theorem 2.9 in [23]). Then due to 
Proposition 3.2 in Magal and Zhao |23j or Theorem 2.3 of Thieme |31j . Lemma 
13.61 ensures that Lemma 13.31 {Hi) holds true. This completes the proof of 
Lemma 13.31 


3.2. Proof of Lemma 13.41 In order to prove that Lemma 13.41 holds true we 
investigate some spectral properties of the Hille-Yosida operator A+DFq {x%) : 
D{A) <Z X ^ X. Note that one has: 


DF^ ( 4 ) 


( 


/ —piiSE{-)oiI — fdjlEP \ 

0 


0 

ai 


ai(3ip*j [sb] -I- PiIeP*i [p] 

\aj) 


\aif3ip*j 4 ] -f PiIePj [p]J 


Hence it is easy to check that (see |5|) 


tUQ.ess DFq {^%)) ^ ~ [m + PiIe] ^ 


where a;o,ess denotes the essential growth rate while the index Xq corresponds 
to the part in Xq of the linear operator. As a consequence of this remark, if 
a {A + DFq {x%)) denotes the spectrum of the linear operator and recalling 
the definition of H in Assumption 13.11 then one obtains using the results in 

[TolIM] that 

fT (A -I- DFq {x%)) n H is only composed of point spectrum. 


As a consequence of the spectral mapping theorem, since the linear semigroup 
{L{t)}t>o is generated by the linear operator {A + DFq {x%))^, the part of 
A -f DFq {x^^ in Xq (see definition in (13.301) 1. in order to prove Lemma [3^ it 
is sufficient to prove that 


(3.37) cr (A -t DFq (4)) n H c {A e H : Re (A) < 0}. 


Let us notice that, since DFq (x^) Xq C Xq, one has {^A + DFq = 

Ao + DFq (Xg)|^^. 

Now it is easy to check that 0 ^ cr (A -|- DFq (x^)). Next let A G H \ {0} and 
X = {ip, 0, S/, oj)'^ G Xq be given. Then consider the problem 

(A - (Axo + DFq (xe)))) {(f,0,a/,aj)'^ = x. 


{ A(fi + (fi' + {/z + (3iIe)'^ + l3iaE{-)ai = ip, 

m = 0 , ^ 

Xai = PjIe pi{a)ip{a)da + a}, 

Xaj = PjIe pj{a)ip{a)da + a}, 


It re-writes as 
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(ip{a) = Jg e (^+M+/3j‘-fE)(“ ['tjj(a') —/3iSE{a')ai]da', 

< A(A)a/ = PjIe pi{a) /“ + op, 

[Xaj = PiIe pj{a)(p{a)da + a} 


As a consequence of the above computations, if A € satisfies A(A) ^ 0 then 
X G p (A + DFq ). 

On the other hand, let A G O be given such that A(A) = 0. Since since A 7 ^: 0 
then one can easily check that 


0 

1 

\aj,x) 


G Ker (a — (A + DFq {xe))q) 


wherein we have set 


' Pi 

(fx{a) := -^je“^“sE(a) / da' and aj^x = ^Ie / pj{a)(fx{a)da 

Jo Jo 

As a consequence, one deduces that 


cr (A + DFo {x°e)) n O = {A G O : A(A) = 0}, 

and (13.371) follows from Assumption 13.11 {in). This completes the proof of 
Lemma 13.41 


3.3. Proof of Lemma 13.51 The proof of the estimate stated in Lemma 15751 
will follow from Gronwall’s inequality. The proof of (13.281) requires some ad¬ 
ditional notation. For each k > 0 we set M{k) = ^ + k. Recalling Definition 
(13.271) and B = Bxo (0; A7(k)) D Ao+, we introduce the quantities: 

Nfo{^) ■■= llAllup. B and Ng{^) = ||G||Lip, 5 ■ 

As a first step let us prove the following lemma: 

Lemma 3.7. There exist some constants M** > 0 and p > 0 such that for 
each K > 0, each e > 0 and each {x, y) G (0) G ^ 0 +) ■ 

(3.38) \\US)x-US)y\\ < Vt > 0 . 

Proof. Let k > 0 and e > 0 be given. Then for each {x, y) G {Bxq (0, k) D Ao+)^ 
one has: 

(3.39) 

Ue{t)x - Ue{t)y = TxoiOix - y) 

-I- lim [ TAo{t-s)XR{X-,A)[Fs{Ue{s)x) - Fs{Ue{s)y)]ds. 

>^^00 Jo 

Now recalling estimate {ii) in Theorem 12.31 one gets that for each s > 0: 

Ue{s)x, Ue{s)y G Bxoi0,M{K))r\Xo+. 

Hence we infer from (13.391) that for each t > 0: 

\\UXt)x - Ue{t)y\\ < ||rAo(t)||||a:-y|| 

limsup A ||i?(A; A)||£(^) [Nfo{k) + eNoiK)] f \\TAo{t - s)||||t/e(s)a: - 17e(s)y||ds. 
A—voo Jo 
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Recalling that operator A satisfies the Hille-Yosida property, the estimate 13.381 
stated in Lemma 1X71 follows from Gronwall’s inequality. q.e.d. 

Now note that since Fg is of class one gets that: 

(3.40) £m ||i7o|lLip,s,„(,o^,,)nx„+ = « with Hg := Fg(.) - DFg {x%) . 

Now let us fix tg > 0. Recalling that Ue{t)x‘% = one infers from Lemma 
EH that there exists a constant M** > 0 such that for each 5 € (0,1), for each 
X S Bxo I”' "^^0+ and each e G [0,1]: 

\\Ue{t)x -x%\\ < - XeW < Vt G [0,fo]- 

As a consequence, since x^^ ^ x^ as e ^ 0, there exists <5o > 0 small enough 

such that for each 6 G (0, i5o) and all e G [0, <5] and all x G Bxq S) fl Xg+: 

(3.41) \\Ue{t)x - x%\\ < M^^:S + sup ||x|;-x^||, VtG[0,to]- 

eG[0,i5] 

Next let S G (0, Sg) be given. Let x,y € Bxo G ^" 0 + be given. Then one 

has for any z G {x, y} 

Ue{tg)z-L{tg)z = lim [ TAo{t-s)XR{X; A) [fE {Ue{s)z) - DFg (x^) L{s)z] ds. 

A^-OO Jg 

If we set 

A®(t, X, y) = [Us{t)x - L{t)x] - [Ue{t)y - L{t)y ], 
recalling the definition Hg in (I3.40|) one obtains 

A^(t,x, y) = e lim [ Txait - s)XR{X; A) [G {Ue{s)x) - G {Ue{s)y)] ds 
>'->-^Jg 

+ lim [ Tao (t - s)XR{X; A) [Hg {Ue{s)x) - Hg {Ue{s)y)] ds 
Jg 

— lim [ TAo{t — s)XR{X;A)DFg (x%) A'^{s,x,y)ds. 

X^OoJg 

If we set K{S) = + supEgfo^^j \\x% — x^||, one obtains due to Lemma EH 

and (13.411) that there exists some constant C > 0 (depending on tg and (5o) 
such that for each t G [0,to]: each e G [0,(5]: 

||A®(t,x,y)|| < Cllx-yll [e + m{S)] + G [ ||A^(s, x, y)||ds, 

Jo 

wherein we have set m(,5) := ||i7o(-)llLip,:Bxo(x«;.if(5))nXo+- Hence Gronwall’s 
inequality applies and provides that for some constant G: 

l|A^(to,a:,y)|| < C[e + m{6)] ||x-y||, 

for all x,y G Bxo (xet^) G Ao+, for all d G (0, i5o] and for all e G [0,(5]. As a 
consequence, one obtains that for all S small enough 

sup \\Ue{tg) - L{tg)\\ j^^^„ ^ < G [(5 + m((5)] . 

£e[0.<5] \ E’ J u+ 

Recalling that Ar(^) —>■ 0 as 5 —>■ 0+ and (13.401) . we conclude that m{6) —>■ 0 as 
(5 —> 0+ and the result follows. 
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4. Application to the prototypical function (11.81) 


In this section we will check that Assumption 13.11 for the prototypical func¬ 
tion (US is indeed satisfied. We will check (ii) before going to {in). 


4.1. Assumption [3711 (ii). In this subsection we will derive some conditions 
on the parameters in order for Assumption 13.11 {ii) to holds true with the 
prototypical function (11.81) . To do this we first give the following general result: 


Proposition 4.1. In addition to i?o[0] > 1; "we assume that function pi S 
T“(0, oo) (d oo) satisfies 

(4.42) p'j{a) < {p + l3ilE)pi{a) a.e.. 

Then Assumvtion lS. 1\ (ii) holds true. 


Before proving this proposition, let us re-write the hypothesis (I4.42p with 
the prototypical function (11.81) . With such a function, pj takes the form 

(4.43) Pi{a) = 1 — Ke~^°‘, Va > 0. 


for some k G (0,1) and r > 0 . Together with this function, note that Rq = 
i?o[0] reads as 


Rn — A 


pLVl 


1 — K- 




p + r 


When i?o > 1, recall that A^; is uniquely defined by the resolution of the 
equation 


1 = A 




vi{p + Xe) 


1 — K- 


^ -l- A 


E 


Xe + tL + r 


Hence ()4.42l) re-writes as the following inequality 


[r + {p, + A_e)] kc ™ < (/r -I- A_e), Va > 0. 


It is then enough to satisfy 

(4.44) ^<(^ + Ai5). 

This condition re-writes as 

A—p*j < 1 < i?o- 

VI 

The former condition means that using function (11.81) . Assumption 13.1| iii) is 
satisfied when the endemic value A^ is sufficiently large. 

It now remains to prove Proposition l4.ll 


Proof of Proposition \4-l\ The proof of this result follows some ideas developed 
by Magal et al in |19j . In this proof, for notational simplicity, we shall denote 

x% = (s_E, 0, Ie, Je) ■ 

Now let us first notice that coupling bounded dissipavity and asymptotic 
smoothness in Theorem [53] together with uniform persistence for Uq in Lemma 
13.31 the results in |23| ensures the existence of an interior global attractor 
for Uo, denoted by Aoo C Mq (see (13.251) for the definition of Mq). 

Now we shall construct a suitable Lyapunov functional on Aqo in order to 
prove that Aqo = To that aim let us consider the function g : x x — 
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lii(a;) —1. Let x G ^oo be given and let us consider {{s{t ,.), J(t))^}tgR C 

^00 an entire solution of the semiflow {l/o}{i>o} ■); 0; '^(0))^ 

X. Next let us consider the uniformly bounded map 

(4.45) V[s,0,I,J]{t) = J a{a)g da + lEg , 


wherein we have set 


a{a) = pi{a)sE{a). 
Now let us notice that one has 


(4.46) 


djS r se' 
S£; L s - 



+ /?/ 


S - Se 
Se 


[Ie 


/], Vt e R. 


Next this yields the following: for each t £ R: 

(4.47) ^ [s, 0 ,/, J] (t) = - y a{a)dag da = J a'{a)g da. 

Now note that (14.421) re-writes as a'(a) < 0 for almost every a > 0. To complete 
the proof of Proposition 14.11 let us consider an increasing sequence {tn}n>o 
such that tn —>■ —oo as n —>• oo. Consider also the sequence of shifted maps 
(s„, 0, In, Jn) {t) ■= {s, 0, 1, J) {t + tn). Then possibly along a subsequence, one 
may assume that 

(.Sn, b; I-n, In') (^) ^ ^S,0, I, (t) aS Tl ^ OO, 

locally uniformly with respect to t £ R with value in Xq-i-- Moreover due to 
(I4.47|) and a' < 0, one concludes that 

[ a'{a)g f ^ da = 0, Vt £ R. 

Jo \SE[a)J 

We infer from the former equality that 


s(t,.) = SE, Vt £ R. 

Furthermore function ^/, satisfies: 

f (da + g) SE^a) = -l3iI{t)sE{a), t £ R, a > 0, 
S = Ht) iPiP*! [se(-)] - ^i) = 0, 

[ J'(t) = PiT{t)p*j [se(-)] - t G R- 

As a consequence, one obtains that 

(/, j) (t) = {Ie, Je) , Vt £ R. 

In addition, this yields 

lim V [s, 0, /, J] (t) = V [sB, 0, Ie, JeJ = 0. 


Finally, since V [s, 0, 1, J] {t) > 0 and since this map t ^ C [s, 0, J, J] {t) is 
decreasing, one concludes that 

V[s,0,I,J] {t) = 0, Vt £ R, 

so that {s,0, 1, J) {t) = {se,0,Ie,Je) for all t £ R and x = x%. This shows 
that Aqo C {x%} and completes the proof of Proposition l4.ll q.e.d. 
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4.2. Assumption 13.11 {iii). In this subsection we shall discuss Assumption 
Id.ll fiiil with a prototype function (11.81) . namely p/ : a i—>■ 1 — . In such a 

case function A : C (see Assumption l3.1D reads as 

/■oo 1 _ —\a 

A(A) = A + pjlE / ^-da, VA e II. 

Jo A 

Then our next lemma reads as: 


Lemma 4.2. For each k S [0,1] the equation 

A G n and A(A) = 0, 

only has roots with strictly negative real parts. 

Remark 4.3. One can see A(A) by extension as: 

poo pa 

A(A) = A + /3|/b / / e-^^’da'da, VA G H. 

Jo Jo 

Proof. Let us first notice that there exists p > 0 and M > 0 such that for 
all K G [0,1]: 

(A G C and A(A) = 0) ^ (p < |A| < M). 

Now note that for k = 0 the equation only has roots with strictly negative real 
parts. Indeed for k = 0 the equation reads as 

1 1 
/r + Xe p + a + Xe 

Hence recalling that Xe = IePi, the equation A(A) = 0 reduces to 
if- + Xe) A^ + (p + Xe)^ X + PjKXe = 0. 


A(A)=A + ^^ 


Thus for K = 0 the characteristic equation only has roots with strictly negative 
real parts. 

To complete the proof of the lemma we apply Rouche’s theorem. It is 
therefore sufficient to show that for each k G [0,1] the equation A(A) = 0 does 
not have any root on the imaginary axis. 

Let K G [0,1] be given. Let us argue by contradiction by assuming that there 
exists w > 0 such that A{iuj) = 0. Then note that this re-writes as: 


—-I- lE0jA 


1 


fi Xe M T i:^ TA^; r p. Xe 
Taking the imaginary part and since w > 0 yields: 

1 K 


r + p + iu: + Xe 


= 0 . 


{p + Xe)^+u}^ {r + p + Xe)^ + 


= 0 . 


Now the above equality is impossible since k G [0,1] and r > 0. This completes 
the proof of the lemma. q.e.d. 
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5. Concluding remarks 

In this work, from two systems with age of infection and/or chronological 
age, we derived an age structured epidemic system with additional mortalities 
and aggregated infectives modelling the dynamics of transmission of a disease 
like Hepatitis B virus. The classical forward bifurcation is retrieved: 

(i) If i?o < 1, then the system with aggregated infectives has a unique global 
asymptotically stable disease free equilibrium. Then the disease dies out. 

(ii) If i?o > 1 and we assume that chronic carriers (most of time asymp¬ 
tomatic) have a very low infectious rate, the system with aggregated 
infectives has (with strong uniform persistence) an unstable disease free 
equilibrium and a global asymptotically stable endemic equilibrium. 

The proof of (ii) uses the useful and complex Theorem 1.2 derived by Magal 
in m- it shows that it is intuitive but technically difficult to claim that the 
global stability property of a steady state remains after a slightly perturbation 
of a parameter. 

The result (ii) was an open problem presented in the discrete age formulation 
of the model studied herejU Fig. 3, page 62]. Consider a continuous age model 
provides a more realistic study of the long term behaviour of the hepatitis 
B disease [38] . Moreover the uniform strong persistence of the semiflow is 
obtained. Roughly speaking here for i?o > Ij uniform persistence is the notion 
saying that the closed subset of extinction for the populations of (acute and 
chronic) infectives is repelling for the dynamics on the complementary set |20| . 
In some sense the persistence points out then the long term survival of infectives 
over sufficiently large time. 
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